Methods for Producing Optical Coherent State Superpositions 
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We discuss several methods to produce superpositions of optical coherent states (also known as 
"cat states"). Cat states have remarkable properties that could allow them to be powerful tools for 
quantum information processing and metrology. A number of proposals for how one can produce 
cat states have appeared in the literature in recent years. We describe these proposals and present 
new simulation and analysis of them incorporating practical issues such as photon loss, detector 
inefficiency, and limited strength of nonlinear interactions. We also examine how each would perform 
in a realistic experiment. 
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I. INTRODUCTION 



For many years a primary goal for researchers in quan- 
tum optics has been the generation of exotic quantum 
states of light. These have included squeezed vacuum 
states, photon number eigenstates (also known as Fock 
states), and many examples of entangled sets of pho- 
tons. States such as these have served well in many ba- 
sic tests of the foundations of quantum theory, and they 
may eventually prove to be useful for quantum computa- 
tion, quantum communication, quantum metrology, and 
lithography. [l[ Therefore there is much interest in im- 
proving methods to produce these states and to generate 
other types of optical quantum states. 

The particular quantum state of interest in this pa- 
per is a superposition of two coherent states with oppo- 
site phase, which is often referred to as a (Schrodinger) 
cat state. Cat states may be used as the logical qubit 
basis in a quantum computer 0, ||. They may also 
serve as input states to an interferometer that is able to 
measure distances with greater accuracy than achievable 
within the limits usually imposed by the light's wave- 
length 4]. Transforming a single coherent state into a 
cat state through unitary evolution alone would require 
a strong nonlinearity. Also, cat states are extremely 
sensitive to decoherence from photon absorption. For 
these reasons cat states containing more than one photon 
on average have been produced only in cavity-quantum- 
electrodynamic experiments in which an atom interacts 
with the electromagnetic field confined to a high finesse 
optical cavity d, Q . In experiments of this sort the cav- 
ity confines the optical mode to a small volume so that it 
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interacts very strongly with an atom passing through the 
cavity. Unfortunately, because the cat state is confined 
to a cavity, it can neither be manipulated with tools such 
as beam splitters or phase shifters, nor be measured with 
standard optical means such as photon counters or ho- 
modyne detection. For the uses described in 0, [H, 0] we 
require cat states that occupy freely propagating optical 
modes. 

In recent years researchers have proposed several 
schemes to produce freely propagating cat states. The 
purpose of this paper is to provide a comprehensive re- 
view and critique of these schemes. We examine the 
performance of these schemes in realistic experimental 
environments subject to problems such as photon loss, 
detector inefficiency and noise, and limited strength of 
nonlinear interactions. We also make recommendations 
for the design of cat production experiments. In this in- 
troductory section we review the properties of cat states 
and optical tools such as beam splitters and photon coun- 
ters, which are commonly used in the cat production 
schemes. The introduction also describes our model for 
photon absorption, how photon absorption affects cat 
states, and the difficulties of verifying that a cat state 
has been produced in an experiment. Section 2 exam- 
ines a method originally proposed by Yurke and Stoler 
0] to transform a coherent state into a cat state using 
the optical Kerr effect. Section 3 briefly discusses the 
suggestion by Wolinsky and Carmichael [8| that one may 
make cat states using a degenerate parametric oscillator 
in the strong coupling regime. Section 4 examines the 
method proposed by Song, Caves and Yurke [9( to pro- 
duce cat states based on an optical backaction evasion 
measurement. Section 5 examines the method of photon 
subtraction from a squeezed vacuum state proposed by 
Dakna and co-authors (To| . Section 6 examines a method 
from Lund, Jeong, Ralph, and Kim [ll| by which one may 
use small amplitude cat states, linear optical devices, and 
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measurements to produce large amplitude cats. Section 7 
examines two schemes [lU, Il3j that use a weak Kerr effect 
followed by measurements to make high amplitude cats. 
We conclude with some general observations in Section 
8. 

A few experiments have already demonstrated the 
production of small cat states containing less than one 
photon [IJ, [TH, [l?]]. These are sometimes called 
"Schrodinger kittens" or simply "non-Gaussian states" . 
All of these experiments use the photon subtraction 
method and will be discussed in Section 4. 



A. Properties of Cats 

For a general introduction to the quantum properties 
of light, we refer the reader to the text [l8|. The elec- 
tromagnetic field may be decomposed into independent 
modes, each of which is defined by a particular polar- 
ization, distribution of light frequencies, and direction of 
propagation. Each mode is modeled as a quantum me- 
chanical oscillator whose frequency is equal to the light's 
frequency u>. The number of photons occupying a mode 
is equal to the number of quantized energy excitations in 
that mode's oscillator. Let a be the photon annihilation 
operator, so a) is the creation operator, and a^a = n 
is the operator corresponding to the photon number ob- 
servable. We use \n) as the photon number eigenstate, 
so n\n) = n\n). For cases in which we are interested in 
more than one mode, we will distinguish the operators 
and states with subscripts labeling each mode. 

The coherent state \a) is the eigenstate of the annihi- 
lation operator with eigenvalue a: 

a\a) = a | or), (1) 

where a may be any complex number. The light pro- 
duced by ordinary lasers well approximates a coherent 
state. (In fact it is more accurate to say that laser light 
approximates the mixed state J d<\>\e % ^ 'a) (e 1 ^ 'a\, but if all 
light modes are phase locked to a local oscillator, then 
this distinction is irrelevant.) In the photon number basis 
a coherent state has the decomposition 



a 



-M 2 /2 



2^ ~7T\ 



(2) 



For ease of notation we usually assume a is a real, pos- 
itive number. In a few cases we will multiply a by a 
complex phase to characterize coherent states with com- 
plex amplitudes. 

We define a cat state to be a superposition of two co- 
herent states with opposite phases: 



|*±(a)) = 



1 



(\-a)±\a)). 



(3) 



The normalization factor N±(a) — 2±2e 2a is required 
because the states |— a) and \a) are not orthogonal. How- 

so that for a moderately sized 



a = 2 this overlap is only 1.1 x 10 1 . When expressed in 
the photon number basis ^±(0)) becomes 



-a 2 /2 °° n 

|*±(«)) = ^=== £ -= ((-1)" ± 1) 



(4) 



| *+(«)) contains only even numbers of photons, and 
I*- (a)) contains only odd numbers. For this reason they 
are often called the "even cat" and "odd cat" states. The 
even and odd cat states are orthogonal to one another, 
and one may perform a measurement to distinguish the 
two by counting the photons. 

In addition to the photon number basis we may also 
measure the cat state in the quadrature bases correspond- 
ing to the position and momentum of the oscillator. The 
position observable is defined by 



x = —= (a 1 + a) . 

V2 V ; ' 

and the momentum observable is 



■1 



(5) 



(6) 



In the x basis, the cat states |^±(a)) have the wave func- 
tions 

(7) 

It has two characteristic Gaussian shaped humps located 
at x = ±\/2<x In the p basis the cat states have the wave 
functions 



*£>(p) 



27T- 1 / 4 



e- p I 2 cos(V2pa) 



e p / 2 isin(V2pa). 



(8) 
(9) 



These are characterized by oscillations of frequency \/2a 
inside a Gaussian envelope centered at p = 0. Resolving 
these oscillations is another possible method to distin- 
guish between even and odd cats. 

We will also consider slightly more general cats with 
phase 0: 



1 



J4>\ 



(10) 



-a\a) 



-4c* 



|*o(o:)) = |#+(a)) and I **■(<*)) = !*-(«))• Changing cj> 
changes the number of photons contained in \^^(a)) and 
requires some type of nonlinear interaction. 



B. Beam Splitters 

Many of the following cat state production schemes 
use beam splitters, so we will introduce them here. For 
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a more thorough discussion of beam splitters see [19j |. 
One may think of a beam splitter as a partly silvered 
mirror. A mode of light enters the beam splitter on each 
of the mirror's two faces. Some fraction T of the light 
energy in each mode is transmitted through the mirror 
and R = 1 — T is reflected. While passing through the 
beam splitter, the state of the two input modes will evolve 
according to the transformation 

MiMa-^CTMiMa, (11) 

where B(T) is the unitary operator given by 

B(T) — e arccos (v / T)(aia2 _a l a 2) 0-2) 

When two coherent states \a)i and | /?) 2 enter a beam 
splitter, their state will become 

|a)i|/3) 2 ^ av / T-/Vl Z ~T) pVT + aVT^ . 

It should be clear from this relation that using only co- 
herent states and beam splitters one can never produce a 
cat state, because a beam splitter simply transforms co- 
herent states into coherent states of different amplitudes. 



C. Photon Counters 

To make cats we will imagine that we have devices that 
can count the number of photons in an optical mode. 
Such devices have been very difficult to engineer, but 
significant progress has been made in this area during re- 
cent years. Photon counters suffer from inefficiency and 
dark counts. Inefficient detectors fail to register all pho- 
tons arriving in the mode of interest. We model this 
by imagining that the mode of interest passes through a 
beam splitter whose transmissivity is equal to the detec- 
tor's efficiency r\ before entering an ideal detector [ijj]. If 
the signal mode has a probability distribution to measure 
n photons given by P(n), then the detector will have a 
probability of P n (jn) to register m photons: 

PuH - £ hf(l - v) n - m P(n), (14) 

when m < n, and P v (m) — when m > n. 

In addition to inefficiencies, detectors also suffer from 
dark counts. These are events in which the detector reg- 
isters more photons than in fact arrive in the mode we 
are measuring. To model the dark counts we imagine 
that a second light mode (in addition to the mode we are 
measuring) enters the detector. We assume the number 
of photons in this second mode has a Poisson probabil- 
ity distribution. If the mean number of dark counts is 
d during a detection event of fixed duration, then the 
probability that q dark count photons are registered is 



Given an initial photon distribution in the signal mode of 
P(n), we combine the dark count and inefficiency influ- 
ences to find that the probability to register m photons 
is 



Pd(m) = 2J P v (x)pd(m - x) 

m 



(16) 



(17) 



Pd(q) 



<Pe~ 



Here the sum adds the probabilities that x photons are 
transmitted through the beam splitter and that m — x 
additional dark counts occur for all x < m. 

Photon detectors are usually characterized by their 
dark count rate D, the number of dark counts registered 
per unit time. The size of d will depend on each detec- 
tor's dark count rate and the duration of the observation. 
Shorter observing times will give a smaller d but require 
faster detectors. Our model for photon counters does 
not incorporate the effects of dead time. We assume that 
such effects can be made to be negligible by increasing 
the time between detection events. 

Let us briefly describe two photon counters that may 
be suitable to assist in a cat production experiment: 
the visible light photon counter (VLPC) and the su- 
perconducting transition edge sensitive photon counter 
(TESPC). The VLPC is discussed in [Sj. It is essentially 
an array of avalanche photodiode detectors (APDs), 
which operate by a version of the photoelectric effect in 
which incident light frees an electron from the valence 
band lifting it into the conducting band. A large voltage 
then accelerates this electron, which collides with oth- 
ers, creating an avalanche of electrons that is amplified 
and manipulated with electronic devices. Most avalanche 
photodiodes are unable to distinguish between 1, 2, 3, 
... photons. Because a single photon depletes the avail- 
able electrons, the arrival of a second photon is unable to 
create a larger electronic signal. The VLPCs can distin- 
guish between different numbers of photons because the 
avalanche is confined to a small portion of the detector's 
area, leaving most of the VLPCs area available for the 
detection of the other photons. References 
discuss the use of the VLPC to distinguish between m— 1 
and m photons with m as large as 10 photons. They also 
successfully demonstrate the use of the transformation 
Eq. (TTJ| to correct for detector inefficiencies. In |2CM ] 
Takeuchi and co-authors report achieving a maximum 
single photon quantum efficiency of r\ = 0.882 ± 0.05 us- 
ing 694 nm light, with a dark count rate of D = 2 x 10 4 
counts/second. Lower dark count rates may be achieved 
by reducing bias voltage applied to the photodiodes, but 
this also reduces the efficiency. They also claim that their 
photon counting system roughly obeys the relation 



(15) 



log 10 (D) = 3.22677+1.206, 



(18) 
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where D is measured in counts per second. In [22|, 
Waks and co-authors describe an experiment in which 
the VLPC was used to observe the photon statistics of 
a squeezed vacuum state. They produced squeezed light 
of 532 nm pulsed at 20 ns. By limiting the detection in- 
tegration window to 20 ns they were able to reduce the 
average number of dark counts per detection event to 
d~4x 10" 4 . 

The TESPC detects photons through calorimetry [25[. 
The light is absorbed by the superconducting electrons in 
a small sample of tungsten. The tungsten sample is volt- 
age biased on the edge of the superconducting-to-normal 
transition, so a small increase in the electrons' tempera- 
ture creates a large change in the tungsten's resistance. 
Measuring this resistance change allows one to calculate 
the energy deposited by the incoming photons. Rosen- 
berg and co-authors [25| report an efficiency of r\ = 0.88 
at a wavelength of 1550 nm. They observed a dark count 
rate of only D = 400 counts/second. This is caused pri- 
marily by blackbody radiation entering the optical fiber 
delivering the photons to the detector and ambient light 
scattering into the fiber. Therefore the dark count rate 
should be significantly reduced for light at shorter wave- 
lengths. Because these primary dark count sources are 
not intrinsic to the TESPC, it is likely that the dark 
count rate can be reduced to negligible levels with suffi- 
cient technical care (for example, shielding and filtering) 
of the light path. The recent advent of such high effi- 
ciency, low dark count detectors such as the VLPC and 
the TESPC makes several of the following cat production 
methods just now feasible. 



D. Cat State Decoherence. 

In addition to the difficulty of creating cat states we 
must also face the extreme fragility of these states to de- 
coherence. We expect that the primary source of decoher- 
ence in an optical experiment will be photon absorption. 
This can be modeled by assuming some of the cat state is 
lost via a beam splitter type of interaction. The cat state 
enters one mode of the beam splitter, and the vacuum 
enters the other mode. After passing through the beam 
splitter some of the cat's energy will be transferred to the 
second mode and lost to the environment. Suppose the 
cat state passes through some medium whose transmis- 
sivity is 77, the cat occupies mode 1, and mode 2 contains 
the vacuum used to model the environment. The initial 
state of the cat /environment system is |^ r ±(cc))i|0)2- Af- 
ter the two modes pass through the beam splitter their 
state becomes 



I - a Vv)i\ - a V l - Vh + \ay/rj)i\a\/l - vh, (19) 

where we have omitted the normalization. The second 
mode is lost to the environment, and the first mode is in 




FIG. 1: Plot of F+(a,j]), the fidelity of a cat state |*+(a)> 
that has suffered from some decoherence by passing through 
a medium of transmissivity r\. The plot includes curves for 
even cats with a = 1 (solid curve), a — 2 (small dashes), 
a = 4 (medium dashes), a = 10 (large dashes). 



a mixed state given by 

p±(a,v) = (1-P±)|*±(a0?))<*±(a0?)| 

+P±|*= F [(a v ^><MaV'7)l> (2°) 

where 



1 N^ay/rf) 
± 2 N±(a) 



-2a (1—rj) 



)■ 



(21) 



and N± is defined just following Eq. Absorption 
affects the cat states in two ways. The amplitude of 
the coherent states is decreased from a to a^/rj. With 
a probability P+ the even cat becomes an odd cat, and 
with probability P_ the odd cat becomes an even cat. 
These probabilities increase quickly with increasing a. 
This result is discussed in more detail in (2(| [2?| • 

To compare the decohered cat state with the original 
cat state, we will use the fidelity, calculated according to 



F±{a,n) = ($ ± (a)\p±(a,r))\$±(a)). 



(22) 



We find that 



F±(a,ry) 



fl-P+ 



4 e -a 2 (l+r,) 



(23) 

In Fig. [1] we investigate the dependence of this fidelity 
on a and rj. This plot shows a very large loss of fidelity 
for higher amplitude cat states even for r] very close to 
1. We also see the distinct behavior of the two types of 
error. For large a the fidelity quickly drops to 1/2 as rj 
decreases from 1, because the initial pure cat becomes the 
mixed state l/2|* + (a))(* + (a)| + l/2|#_(a))(#_(a)|. 
As 77 continues to decrease, the coherent state amplitude 
decreases, lowering the fidelity from 1/2 to 0. 
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E. Cat State Verification. 

Any cat generation experiment requires a verification 
that the experiment makes cat states. This is likely to be 
accomplished by performing a series of homodyne mea- 
surements, which can be used to reconstruct the state 
through quantum tomography [H, [28| . If the photode- 
tectors used in the homodyne measurement have effi- 
ciency 77, then the homodyne measurement signal will be 
equivalent to that of the generated state having passed 
through a beam splitter with transmissivity 77. If 77 
is known, the homodyne detection inefficiency effects 
can be corrected in data analysis by using the inverse 
Bernoulli transformation (based on Eq. (fH)) ) and other 
techniques described in [lj| l28j ]. The correction is ap- 
proximate and depends on the accuracy with which the 
efficiency is known, what other error sources are present, 
and the algorithm used to infer the quantum state from 
the homodyne data. Fortunately, the photodetectors 
used for homodyne measurements have much higher ef- 
ficiencies than those used for counting small numbers 
of photons. For example Polzik, Carri, and Kimble re- 
port obtaining an overall homodyne detection efficiency 
of 7? = 0.98 in flU. 

Reconstructing a cat state through tomography re- 
quires a large number of homodyne measurements. In 
this procedure one measures the probability distribution 
of the x quadrature for each phase of the local oscillator 
from to 7T. In practice this is done by constructing a 
histogram of x with a finite bin size, with a finite number 
of measurements, for a finite number of local oscillator 
phases. We would like to estimate the required number 
of phases, the bin size resolution, and the number of mea- 
surements required for each histogram. Determining the 
best method for assigning uncertainty to, for example, 
the density matrix elements of the reconstructed state 
and the relationship between error bars and the number 
of measurements requires further research. For discussion 
of these issues see [2^, [3(J HH . A full treatment of this 
difficult data analysis problem requires further research. 
Instead, we quote a few rules of thumb from [19| . which 
contains a much more detailed discussion. To reconstruct 
the state of a mode containing at most M photons, one 
should collect histograms for M + 1 equally spaced local 
oscillator phases, and the histogram bin size resolution 
must be less than 7r/(2\/2M + 1). Suppose the state we 
are measuring has the photon number probability distri- 
bution P(m). To accurately reconstruct the probability 
to detect m photons each histogram should contain on 
the order of 4 x (P(m)) -2 counts. These rules can help 
us to estimate the order of magnitude of the number of in- 
dividual homodyne measurements required to verify the 
production of a cat state. For example, the even a — 2 
cat contains a mean number of approximately 4 photons. 
The probability to detect 10 photons in this cat state is 
P(10) ~ 0.01. If we are willing to ignore higher photon 
numbers, reconstructing the Wigner function of this state 
will require on the order of 4x (0.01)~ 2 x (10+1) ~ 4x 10 6 



individual homodyne measurements. For the more mod- 
est goal of verifying an a — 1 cat, we could ignore photon 
numbers larger than 4. In this case P(4) = 0.03, so we 
would need ~ 2 x 10 5 measurements. These rough esti- 
mates highlight the need for any proposed experiment to 
be able to produce a large number of cat states in order 
to perform a full reconstruction of the state generated by 
that experiment. 

It is possible to estimate the fidelity of the state pro- 
duced in an experiment with a perfect cat state without 
fully reconstructing the unknown state. See [32l | for an 
application of such a technique in an ion trap system. 



II. KERR EFFECT 

In 1986, Yurke and Stoler pj considered the evolution 
of a coherent state under the influence of the Kerr ef- 
fect. They showed that, under suitable conditions, the 
coherent state evolves into a cat state. 

The Kerr effect causes a phenomenon known as "self- 
phase modulation." When light passes through a 
medium exhibiting the Kerr effect, the rate at which 
the light's phase advances depends on the intensity of 
the light. This is modeled with an anharmonic-oscillator 
Hamiltonian of the form 

H K = ujh + xn 2 , (24) 

where lj is the light frequency, and \ is the strength of the 
nonlinear term. The first term causes only linear phase 
evolution, so we neglect it in the following treatment. 

If we apply this Hamiltonian to the coherent state \(3) 
for a time , 

|/3> ^|* f (/3)) = -±= (| -/?)+#)). (25) 

Due to the i factor, this cat is slightly different from 
those discussed in other sections of this paper, but the 
i cat exhibits the same interesting and useful features. 
Achieving this transformation requires a large nonlinear 
strength x or a long interaction time. 

The Hamiltonian in Eq. (|24[) is naive and cannot cap- 
ture the full dynamics of the light's interaction with the 
nonlinear medium. We assigned a single frequency to the 
light, and only continuous waves of infinite duration can 
have a single frequency. If a pulse of light with a finite 
duration, with say a Gaussian shape, interacts with the 
nonlinear medium the dynamics will be far more com- 
plicated. Because the pulse will have lower intensity (or 
photon number density) at its leading and trailing edges, 
the Kerr effect will have a greater effect on the center of 
the pulse. Any linear or nonlinear dispersive properties 
of the medium will distort the pulse shape. We are ignor- 
ing all of these effects in the following treatment. Instead, 
our goal is to provide the first analysis of the relationship 
between cat state fidelity, Xi an d absorption. 
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To model loss we assume that the number of photons 
decays in time according to 



(h(t)) = (h(0))e-^, 



(26) 



where 7 describes each material's photon loss rate. The 
evolution of a (possibly mixed) state traveling through a 
Kerr medium with loss obeys the master equation 



dp 
df 



-ix In ,p\ + 
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\ap, a) 



2 l^P^} 



(27) 



This equation has been solved by Milburn and Holmes 
in [33]. They used the antinormally ordered quasiprob- 
ability distribution known as the "Q-function." The Q- 
function allows one to represent p as a function over the 
complex plane, and the value of Q{a) is proportional to 
the probability that p "contains" the coherent state \a). 



Q(a) = ~(a\p\a). 

IT 



(28) 



Note that a may be complex. By projecting both sides 
of the master equation onto (a\ from the left and onto 
I a) from the right, one obtains a differential equation for 
Q(a). Milburn and Holmes solved this using a power 
series and obtained 



Q loss (a) = e^ 2 - 

q,p=0 



(a/3*) 9 {a*(3f 



q\p\ 



exp 



-itj(j> + q) 



4x 



exp 



XiP-q e ~W + 2i(p-q) 
l + 2i{p-q) 



(29) 



after t = Q\ oss is a function of "f/Xi which is an 
unitless quantity. 

We calculate the fidelity of the state described by Qi oss 
with the cat in Eq. (|25|) using the relation 



d 2 aP cat (a)Q ioss (a), (30) 



where P C at is the normally ordered quasiprobability de- 
scribing the state in Eq. ([23)1 . The P function allows 
us to represent any density matrix as a diagonal sum or 
integral over coherent states: 



d 2 \a){a\P{a). 



(31) 



We calculate P C at(a) for Eq. (l25l) u sing the methods ex- 
plained in the third chapter of [3J| and find 



-Peat (a) 



E 



2 1 — ' mm 

m,n— 



(l + (-1)™+™ + ie- 2 ^ 2 [(-1)" - (-1)" 1 ]) 
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da 71 d(a*) 



5 2 (a), (32) 
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FIG. 2: (a) Fidelity of the state produced using the Kerr 
effect in the presence of loss when our goal is to produce 
a cat of amplitude (5 — 1 (solid curve), 2 (long dashes), 3 
(medium dashes), and 4 (short dashes). The amplitude of the 
input coherent state has been optimized to give the maximum 
fidelity, (b) Same as (a), zoomed in on small j/x- 



where S 2 (a) is the Dirac delta function on the complex 
plane. 

We can now evaluate the integral in Eq. (|3T)|) and find 
the fidelity with which one can make cat states using a 
lossy Kerr effect medium. Plots of the fidelity as a func- 
tion of j/x arc in Fig. O Creating a cat of /3 = 1 with 
F > 0.5 requires j/x < 1-5. Larger cats are much more 
sensitive to loss. Making high fidelity cat states directly 
using the Kerr effect requires some nonlinear optical ma- 
terial/process for which 7/x < 1- 

Materials exhibiting the Kerr effect are often described 
using a "nonlinear index of refraction" ri2, which depends 
on the light intensity /. 



n(/) = ni + nil = rt\ + n% 



hujh 
A cS T 



(33) 



where n\ is the usual linear index of refraction, A e g is the 
effective cross-sectional area of a pulse of light traveling 
through the material, and T is the pulse's duration. To 
relate ri2 to x we use 



Uuj 2 n2 



(34) 
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from [35j . To obtain a large 7/x for a given medium, we 
should confine the mode's A e g to be as small as possible. 
One way to achieve this would be to use short pulses 
of T ~ 1 fs traveling through single mode fibers, which 
typically have A e g ~ 7^tm 2 . 

Fused silica fibers used for telecommunication have 
exceptionally low loss of 0.2 dB/km at the wavelength 
1550 nm. This corresponds to 7 = 1.79 x 10 5 s _1 . Fused 
silica has a small n 2 = 2.6 x 10~ 20 m 2 /W [36]. This 
corresponds to x = 620 s _1 for these pulses. Therefore 
single mode fused silica fibers can give 7/x ~ 260, which 
is much too large for making cat states. 

Chalcogenide (AS2S3) glass fibers have larger n 2 — 
2 x 10~ 18 m 2 /W. However, their loss is 100 dB/km 
at 1.55 jtim. Single mode chalcogenide can have 7/x ~ 
1.3 x 10 4 . We have been unable to find any conventional 
nonlinear optical fibers that can give a smaller 7/x than 
that of fused silica [37| • 

Mecozzi and Tombesi [H[ and Yurke and Stoler [3^| 
published papers in 1987 that proposed making cat-like 
states using a two mode cross-phase modulating or four 
wave mixing Kerr effect [H, |39| . They described the 
application of the Hamiltonian with the nonlinear inter- 
action term 

H in t = X (ai&l + a\d 2 ^j (35) 

to the input state |0)i|/?)2- After time 7r/(4x), the two 
modes will evolve to 

- e-^ 4 |0) 1 |*_(/3)) 2 + |vE< + (*/3))i|0) 2 . (36) 

Here the cat may be found in mode 1 or 2 with prob- 
ability 0.5, while the other mode contains the vacuum. 
Tombesi and Mecozzi presented further analysis in [40| . 
where they discuss both self-phase modulating and cross- 
phase modulating Kerr effects. 

A similar scheme was proposed in 1997 by Vitali, 
Tombesi, and Grangier [41j. They described using the 
interaction Hamiltonian 

Hint = X"l"2, (37) 

and the input state |a)i|/3)2 to make 

|¥ + (a))i|/3>2-|*-(a)>i|-/3>a. ( 38 ) 

Then one should measure the ^-quadrature of mode 2. 
With probability approximately 0.5 the measurement re- 
sult will be near ±\/2/3, and mode 1 will be the cat 
l^±( a ))- Vitali, Tombesi, and Grangier also examined 
this scheme's performance in the presence of loss, con- 
cluding that the loss should be very small in order to 
observe the interference fringes in the q-quadrature of 
the cat state. 

We have not performed a full analysis of photon ab- 
sorption in the cross-phase modulating Kerr effect. How- 
ever, the cat states generated by this method are just as 
fragile as those generated with the self-phase modulating 



Kerr effect. Therefore we expect that the requirement 
for a low ratio of loss to nonlinear coupling strength in 
cross-phase modulation is similarly demanding as in the 
self-phase modulation case. 

Two technologies that enhance optical nonlinearities 
are photonic crystals and electromagnetic induced trans- 
parency (EIT). Photonic crystals may be able to reduce 
the group velocity of light pulses by factors of 100 to 
1000 (42l. |43| . A 300-fold reduction of the group velocity 
(compared to c) of 1550 nm light was observed in a silicon 
photonic crystal [44j . Reducing the group velocity would 
allow one to make an equal size cat with a smaller length 
of fiber and a proportional reduction of 7. However pho- 
tonic crystal waveguides such as that used in [44| have 
an attenuation of ~ 3.6 x 10 5 dB/km (45|, which is much 
too large to maintain cat state coherence. 

EIT is a technique for coupling two (or more) light 
fields through an atom [46| . By engineering the inter- 
action between the light fields and the atom, one can 
obtain many interesting linear and nonlinear optical ef- 
fects, one of which is cross phase modulation between 
the light fields [13] • During cross phase modulation, 
the phase of each light mode advances at a rate pro- 
portional to the number of photons in the other mode. 
If two coherent states \a) and \0) are subjected to cross 
phase modulation for the correct time, they will evolve 
to \a, {3) + \-a,/3) + \a, -(3) - \-a, -0) A8\. This is a 
state of two cats entangled with one another. However, 
cross-phase modulation at the level of only a few photons 
has not yet been demonstrated [46] ]. 

III. DEGENERATE OPTICAL PARAMETRIC 
OSCILLATOR 

In 1988 Wolinsky and Carmichael presented an analy- 
sis of degenerate optical parametric oscillators (DOPO) 
8]. In a DOPO a strong pump field of frequency 2w is 
used to pump a nonlinear crystal. The crystal exhibits a 
susceptibility and can oscillate at frequency to when 
driven by the pump field. Through this process pump 
photons are converted into pairs of signal photons with 
frequency ui. To increase the interaction strength the 
crystal is placed in an optical cavity. Photons of both 
frequencies may leak out of the cavity or be absorbed by 
the crystal. 

Wolinsky and Carmichael examined the behavior of a 
DOPO using the Hamiltonian 

H = ^ (a\ 2 a p - a 2 6t) + 1E (fit _ fi p ) + H loss , (39) 

where a s is the annihilation operator of the signal mode, 
and a p is the annihilation operator of the pump mode. 
H\ oss describes the linear losses from the cavity mirrors 
and in the crystal. When the pump is treated as a classi- 
cal field, and loss is neglected, this Hamiltonian produces 
squeezed light. Wolinsky and Carmichael analyzed this 
system, including loss, in the strong coupling regime with 
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g ~ 1 , which is much larger than that achieved in current 
squeezing experiments. Their analysis was based on the 
positive- P representation, an approximation of the usual 
normally ordered P representation. Their analysis sug- 
gested that in the regime where the ratio of the pump 
strength A to g 2 is very small, one may find the cat state 

y^X/g) + — V^/g) in this system. 

However, Reid and Yurke [49( showed that for any fi- 
nite number of photons in the signal mode, the steady 
state of the degenerate parametric oscillator in the strong 
coupling regime is instead an incoherent mixture of 

VA/g) and —\fX/g) ■ The system's Wigner func- 

/ s Is 

tion is always positive, and no interference fringes would 
be visible in homodyne detection. Further analysis by 
Krippner, Munro, and Reid [50] investigated the short- 
time evolution of this system starting with the signal 
mode in the vacuum state. Their calculations indicate 
that the signal mode may exhibit interference fringes in- 
dicative of a cat-like state, but these fringes only survive 
for times much shorter than the signal mode's cavity de- 
cay time. Unfortunately, the coupling strength re- 
quired to produce cats using this method is vastly larger 
than the strength of currently available materials. Also, 
we do not have a good method to confirm the cat's pres- 
ence, because it is trapped in the optical cavity and will 
decohere when it exits the cavity mirrors. 



IV. BACK-ACTION EVASION MEASUREMENT 

Because of the great difficulty in creating cats directly 
using the Kerr effect, researchers have proposed numer- 
ous methods for making cats based on squeezing followed 
by measurement and postselection. The first of these was 
described by Song, Caves, and Yurke in [9(. Their scheme 
is in Fig. [3J It is based on the back-action evasion mea- 
surement described in [5l|, except that the input state 
is the vacuum, and we measure the number of photons 
rather than one of the quadratures. 

We begin by making an entangled, squeezed state of 
modes 1 and 2 through non-degenerate down-conversion. 
This prepares the state 



Si 2 (r)|0i,0 2 ) =e r 



Oi.Oa 



(40) 



The wavefunction describing the x-quadratures for the 
two modes is 



7* 



exp 



^sq (El,2;2) = (a!i,X2|-S , i2(r)|0i,02) = 

I P — 2r / xj ^£2_\ 2 _ I P 2r ( x±_ , X2_ 

2 e \J3 V2j 2 6 \V2 + V2 



(41) 




FIG. 3: Diagram of the scheme to make cats using the back- 
action evasion measurement. Si2(f) is a nondegenerate down 
conversion crystal, which creates squeezed light in modes 1 
and 2 at frequency uj. The dashed line is a beam splitter with 
transmissivity T. The triangle is a photon counter, and Si(s) 
is a degenerate down conversion crystal that squeezes mode 
1. Not shown are a light beam or beams used to pump the 
down conversion processes and mirrors to redirect the light. 



mode 2, obtaining the result m. This projects the state 
onto the m photon eigenstate of mode 2, which is 



4>m (2:2) 



2 /2 H m (x 2 ) 



-x 2 /2 



v/2 



"ml. 



(42) 



where H m (x) is the Hermite polynomial with index 
m. Last, mode 1 is squeezed in the degenerate down- 
conversion 



S\ (s) = exp 



The final state is 

e s/2 



2 



dx 2 (j) m {x 2 ) 



(43) 



x ip s(1 (Vfxie s + Vl - Tx 2 , Vl - Txie s 



Tx 2 ), 
(44) 



where P(m) is the probability to measure m photons. 
We will use ip m and \ip m ) to denote the output of several 
schemes to make cats by post selecting on counting m 
photons, but in each case ip m will represent a different 
state. 

We have performed the first calculations of the fi- 
delity of the state produced from the back-action evasion 
scheme with the perfect cat state wave function (x) 
from Eq. 



F = 



dxi^\x)^ m (xi) 



(45) 



Modes 1 and 2 then meet in a beam splitter with trans- 
missivity T, which performs the transformation xi — > 
VTxi + V 1 — Tx 2 and x 2 — > y/1 — Txi — VTx 2 on 
ip S q(xi, x 2 ). We next count the number of photons in 



To maximize the fidelity with which we can make a par- 
ticular cat, we optimize over the two squeezing parame- 
ters r and s, the beam splitter transmissivity T, and the 
number of photons detected m. [In the original scheme 
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of Song, Caves, and Yurke, T = cos 2 (l/2 arcsin(tanhr)), 
which is the value of T required for a back-action evasion 
measurement, but here we treat T as a free parameter.] 
For a particular to, Wegener ally find a significant region 
of the (r, s, T) parameter space that gives the maximum 
fidelity. 

Table Q] shows the fidelities with which we can make 
an even cat state with a = 2, along with the parameters 
that give the maximum fidelity. Particularly we choose 
to show the combination (r, s, T) that gives the highest 
P(m) with which one can make the highest fidelity cat 
states. By increasing T and reducing r and |s|, one can 
obtain cats with the same fidelity but lower probability. 
The table shows fidelities only for even to, because for 
odd to the fidelity is always 0. 

TABLE I: Fidelities for Making Cats with Back- 
Action Evasion 



rn 


F 


P(m) 


r 


s 


T 


2 


0.9709 


0.110 


1.14 


-1.35 


0.808 


4 


0.9978 


0.056 


1.44 


-1.48 


0.710 


6 


0.9995 


0.038 


1.61 


-1.63 


0.652 


8 


0.9998 


0.029 


1.76 


-1.77 


0.616 



Fidelities with which one can make even cat states with a = 2 
using the back-action evasion scheme. We optimize r, s, and T 
first, to give highest fidelity, and second, to give highest 
probability. 

Let us discuss the case of m — 2 in more detail. The 
values shown in the table show the highest probabil- 
ity with which one can obtain the highest fidelity cat 
states. We can reduce the demand for high squeezing 
by accepting a lower success rate. For example, us- 
ing (r = 0.00235, s = -0.588, T = 0.999999), we ob- 
tain cats with the same F = 0.9709 with probability 
P(2) = 4 x 10~ n . A reasonable compromise could be 
(r = 0.074, s = —0.593, T = 0.999) at P(2) = 4x 10~ 5 . 
This gives flexibility to reduce the required level of 
squeezing while suffering a lower success probability. 

The original back-action evasion scheme of Song, 
Caves, and Yurke was improved by Yurke, Schleich, and 
Walls in [52l |. They recommended squeezing mode 1 be- 
fore the two mode squeezing S\ 2 (r). This removes the 
need to squeeze mode 1 after the beam splitter used for 
the back-action evasion measurement, and it reduces the 
need for strong squeezing. Figure 2] shows the new pro- 
cedure. 

We analyzed this improved scheme using methods sim- 
ilar to those just described for the original back-action 
evasion scheme. We again show a table of fidelities for 
making an even cat state with a = 2 in Table ITT1 Moving 
the squeezing of mode 1 to the beginning of the pro- 
tocol allows us to produce cats with equal fidelity with 
much lower levels of squeezing. For the to = 4, 6, and 8 
we can remove the beam splitter and count the photons 




FIG. 4: Diagram of the improved scheme to make cats using 
the back-action evasion measurement. It is similar to that in 
Fig. [3] except that Si (r) has been moved to the beginning of 
the network. Not shown are a light beam or beams used to 
pump the down conversion processes and mirrors to redirect 
the light. 

in mode 2 just after the entangling operation of S\ 2 (r). 
One can make cats with fidelities equal to those shown in 
the table requiring less squeezing, by accepting a lower 
P(to). For example one may use the combination to — 2, 
r = -0.0009, s = -0.589, and T = 0.547, to make an 
a = 2 even cat state with probability P(m) = 1.01 x 10 -6 
and F = 0.9709. 

TABLE II: Fidelities for Making Cats with the Im- 
proved Back-Action Scheme 



m 


F 


P(m) 


r 


s 


T 


2 


0.9709 


0.110 


-0.263 


-1.36 


0.972 


4 


0.9978 


0.056 


-0.271 


-1.41 


1 


6 


0.9995 


0.0017 


-0.162 


-0.62 


1 


8 


0.9998 


0.000016 


-0.116 


-0.45 


1 



Fidelities with which one can make even cat states with a = 2 
using the improved back-action evasion scheme shown in Fig. [4] 
Wc optimize r, s, and T first, to give highest fidelity, and second, 
to give highest probability. 

We now introduce a new simplification of the "back- 
action evasion" scheme, which is now quite divorced from 
the original back-action evasion ideas of [5l[. Because 
the beam splitter used to interfere modes 1 and 2 before 
the photon counting measurement shown in Fig. 0] is un- 
necessary for to > 4, we remove it altogether. One can 
always turn a two mode squeezing operation into two 
single mode squeezers followed by a beam splitter [53| : 
Su(r) = B 12 {l/2)S 2 (-r)S 1 {r). Therefore the sequence 
Su(r)Si{s) = B 12 {l/2)S 2 (-r)S 1 (s + r), which is just 
independent squeezing of modes 1 and 2 followed by a 
beam splitter. This leads us to consider a scheme that 
begins with squeezed states in both modes 1 and 2, which 
then interfere at a beam splitter with transmissivity T. 
We count the number of photons in mode 2, obtaining 
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TABLE III: Fidelities for Making Cats with Simpli- 
fied Back-Action Scheme 



rn 


F 


P(m) 


r 


s 


T 


2 


0.9709 


0.110 


0.263 


-1.62 


0.665 


4 


0.9978 


0.056 


0.274 


-1.76 


0.497 


6 


0.9995 


0.038 


0.221 


-1.85 


0.398 


8 


0.9998 


0.029 


0.182 


-1.93 


0.332 



Fidelities with which one can make even cat states with a = 2 
using our simplified "back-action evasion" scheme shown in 
Fig. [5] We optimize r, s, and T first, to give highest fidelity, and 
second, to give highest probability. 



FIG. 5: Diagram of our simplified "back-action evasion" 
scheme to make cats. Modes 1 and 2 begin in the vacuum 
state, and are transformed to squeezed states by Si(r) and 
S2(s). The two squeezed states meet at the beam splitter 
with transmissivity T. We then count the number of photons 
in mode 2, obtaining result m. Mode 1 then contains \ipm), 
which should be similar to a cat state. Not shown are a light 
beam or beams used to pump the down conversion processes 
and mirrors to redirect the light. 



cantly degraded. Experimentalists often report results 
of squeezing experiments using the number of decibels of 
the variance of the quadrature noise above or below the 
variance of the vacuum's noise. The number of decibels 
is related to our squeezing parameter s by 



# of dB = -101og 10 e 



1.69s 



(48) 



result m and preparing the output of mode 1 in the state 
\ip m ), which should be similar to a cat. Figure O contains 
a diagram of this scheme. 

After having been squeezed, modes 1 and 2 are in the 
state 



■ exp 



( Xl e r )' 2 (x 2 e 8 ) 2 



The beam splitter causes the the rotation of xi and x 2 
described above. After measuring m photons in mode 2, 
the output state is 

ipm(xi) = / dx 2 <j} m (x 2 )^ sq (VTxi + v 1 - Tx 2 ) 

x Vs q (vl - Txi - Vfx 2 ) , (47) 

where 4> m (x 2 ) is the photon number eigenstate wave func- 
tion from Eq. (|42[) . We show a table with the fidelity 
with which this state approximates an even cat state of 
a = 2 in Table ITTT1 Again in this scheme we can trade off 
squeezing strength for lower P(m) while keeping the same 
fidelity. We could make an even cat state of a = 2 whose 
fidelity is 0.9709 using m = 2, r = -0.590, s = 0.0010, 
and T = 0.9975 with P(2) = 1.3 x 10~ 6 . Here the squeez- 
ing of mode 2 is so small, that we are tempted to remove 
it altogether. We discuss such a scheme in detail in the 
next section. 

Two difficulties for implementing these "back-action 
evasion" schemes are their need for a high efficiency pho- 
ton counter and strong squeezing. If the photon counter 
misses one photon, the scheme makes a state orthogo- 
nal to the desired cat, so the fidelity would be signifi- 



for the a:— quadrature, and the negative of this for the 
p— quadrature. A squeezing parameter of s = 1 corre- 
sponds to ~ — 8.7 dB of quadrature noise power squeezing 
below the shot noise limit. The most impressive squeez- 
ing experiments for continuous wave li ght that the au- 
thors are aware of have achieved —9 dB [5J] and —10 dB 
[HI of squeezing. In a typical pulsed squeezing experi- 
ment one may obtain —3 dB (s = 0.35) of squeezing [56| . 
Fortunately these schemes allow a trade-off in which one 
may use weaker squeezing to make equally high fidelity 
cats with a decreased success probability. 

Here we have made no effort to analyze any impurity 
or loss in the squeezing, which is likely to have signifi- 
cant effects. For discussion of loss we refer the reader to 
a pair of papers by Tombesi and Vitali [13] and Vitali 
and Tombesi [58j . In addition to loss, these papers de- 
scribe a method for detecting the cat state by exploiting 
the back-action evasion interaction. They imagine pro- 
ducing a cat in a cavity that has much lower loss for the 
signal (cat-containing) mode than for the meter mode. 
The meter mode exits the cavity, and we count its pho- 
tons. This prepares the signal mode in a cat state. We 
could then use the back-action evading interaction be- 
tween the signal and meter modes to measure the signal's 
state by performing homodyne detection of the meter's 
quadratures. Using this method requires that the back- 
action evasion conditions be met, so one could not freely 
optimize T for making high fidelity cats. 

It is also well known that hig hly squeezed states are 
very sensitive to phase noise [541 ]. It would be very chal- 
lenging to incorporate the high efficiency photon counter, 
two strong down conversion processes, and stable phase 
control in a single experiment. 
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FIG. 6: Diagram for the generation of cat states by means of a 
conditional photon number measurement on a beam splitter. 
The down conversion Si(r) creates the single mode squeezed 
state in mode 1. It is input into one port of a variable trans- 
missivity T beam splitter with mode 2 containing a vacuum 
state. A definite measurement of m photons on one output 
port of the beam splitter prepares the state \ipm), which is a 
good approximation to a cat state. 



V. PHOTON SUBTRACTION 

A simplification of the back-action evasion scheme was 
proposed by Dakna and co-authors [io| . It too uses 
squeezing followed by photon counting but accomplishes 
similar goals with only one stage of squeezing. Their 
scheme is depicted in Fig.[5]and works as follows. First it 
requires first the preparation of a squeezed state of light. 
This can be accomplished by sending a laser beam at 
twice the frequency of the desired squeezed state through 
a down-converting crystal. This prepares mode 1 in the 
state 



\^sq) 



r{*S-W)*)l2 



|0) 



= (l-^E^l^l I 2 ")' ( 49 ) 



where A = — tanhr. \ip S q)i has only even numbers of 
photons, so in this respect it is already like the even cat 
state. 

To make a cat, we combine the squeezed state and 
a vacuum state on a beam splitter B\2{T) with variable 
transitivity T. On the mode 2 output port from the beam 
splitter we count the number of photons and obtain the 
result m. The conditional state of output mode 1 is then 



1 



oo 



AT 



\n), (50) 



where 



(n + m)!(l + (-l) n+m ) 



(51) 



and M m = J2 n c n,m (if")™ ^ an 0< ld number of pho- 
tons in is registered in the detector, then c n>m = for 
all even n's, and if m is even, then c nyTn = for all odd 
n's. This condition is required by photon conservation; 
because \ip S q)i has only even numbers of photons, the to- 
tal number of photons emerging from the output ports of 
the beam splitter must also be even. The mean photon 
number of \tp m ) is 



1 

Mr 



71=0 



AT\ ,H 



(52) 



Eq. (|50|) can be broken into two cases: the state re- 
sulting from an even m result (which should be similar 
to an even cat state) and the state from an odd m (which 
should be similar to an odd cat). For m even, Eq. ([50]) 
has only even photon numbers and can be written in the 
simplified form 



Vl'n 



1 



E 



(2n + m)! ( 



AT 1 
2 / 



m n=0 



|2n). (53) 



For AT small, this expression can be further approxi- 
mated as 



hU « | ) + Ari-^|2) + 



(54) 



As m increases, so does the population in the |2) (and 
higher) states compared with the m — situation. Thus 
for small AT the mean photon number increases as m 
increases. 

To find how closely the states \^ m ) match true even 
cats, we calculate the fidelity F = \(^> + (a)\ip m )\ 2 be- 
tween the two states. Plots of the fidelity of |>F+(2)) 
versus AT for various m's appear in Fig. One can 
clearly see that the fidelity is improved and the amount 
of required squeezing is lower when more photons are 
detected. 

We now examine how the maximum fidelity (maxi- 
mized over the values of AT) increases with m. In Fig. [8] 
we plot the maximum fidelity for each m measurement 
when attempting to make even cats with a = 1,2, and 3. 
Fig. [9] shows the AT that optimizes the fidelity for each 
m measurement. 

Making a high fidelity cat state of a given a is best 
achieved by conditioning successful cat production on de- 
tecting larger numbers of photons and using smaller AT. 
Unfortunately, large m detection events occur less fre- 
quently as the amount of squeezing is decreased. The 
probability of detecting m photons is 



P(m) 



/ 1-A 2 
' 1 - (AT) 2 

Int[m/2] 

E 



A 2 T(1 - T) 



1 - (AT) 2 



(rn-22)m 2 (2AT) 2r 



(55) 
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FIG. 7: Plot of the fidelity of the state \ip m ) with |*+(2)> 
versus AT for m=0 (solid curve), m=2 (small dashes), m—4 
(medium dashes), m—6 (large dashes). The upper horizontal 
axis shows the number of decibels of squeezing required when 
T = 1 is used. 



1 

0.9 
0.8 

Fidelity 0.7 
0.6 
0.5 
0.4 



2 4 6 



10 12 14 16 18 20 

m 



FIG. 8: Plot of the fidelity of the state |* m ) with |*+(a)> for 
a — 1 (diamonds), 2 (stars), and 3 (squares). Notice that only 
even m are represented, because odd m events give a fidelity 
of zero. For each point we have numerically optimized AT to 
give the maximum fidelity. 
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FIG. 9: Plot of the product AT that maximizes the fidelity 
shown in Fig. [S] 



FIG. 10: Plots of the probability to detect m photons after 
sending a squeezed vacuum state with squeezing parameter 
A through a beam splitter with transmissivity T. Each plot 
shows the probability obtained using the product AT that 
maximizes the fidelity shown in Fig. [9] Plot (a) shows the 
case in which A = 0.99 and the probability is high. Plot (b) 
shows the case in which A = T. Plot (c) shows the case in 
which squeezing is decreased and T = 0.99. 



Although the state \ip m ) and its fidelity with an even 
cat depend only on the product AT, its probability de- 
pends on A and T separately. P(m) can be increased by 
increasing A while decreasing T (while keeping AT fixed), 
but producing highly squeezed light is difficult. Fig. [TU] 
plots the probability to detect m photons for various A 
and T combinations. 
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FIG. 11: Probability to detect m=2 (solid curve), 4 (small 
dashes), and 6 (long dashes) photons as a function of squeez- 
ing (shown with both A and in decibels). In each case T is 
adjusted so that AT gives the best fidelity to produce a cat 
state with a = 2. 



FIG. 12: Fidelity to make a cat state with a — 2 as a function 
of A when using detectors with efficiency r\ = 0.9. The solid 
curve shows the m = 2 case, the small dashed curve shows 
the m — 4 case, and the large dashed curve shows the m — 6 
case. For each combination of A and m, the beam splitter 
transitivity has been adjusted to give the maximum fidelity. 
These maximum fidelities are equal to the maximum fidelities 
achievable perfect detectors are used. 



Let us examine a few illustrative cases. Suppose we 
desire to create an even cat state with a = 2. The 
squeezed state closest to the a = 2 cat has A = 0.883 
and a fidelity of only 0.588. If we detect m = 2 pho- 
tons from the state prepared using AT = 0.613, we will 
create a cat with a fidelity of 0.891. Detecting to = 4 
using AT = 0.469 produces a fidelity of 0.950. Detecting 
to = 6 using AT = 0.380 produces a fidelity of 0.971. 
Notice that higher fidelity is achieved by detecting larger 
numbers of photons and using smaller AT. We would like 
to select A and T to give the highest success probability. 
In Figure QT] we plot P(m) versus A for detecting to = 2, 
4, and 6 photons. The beam splitter's transmissivity is 
adjusted to keep AT equal to the value giving the best 
fidelity for each photon detection case. 

A weakness of this method for creating cat states is 
its need for high efficiency photon counters. If the de- 
tector misses one photon then this scheme will create a 
state that is orthogonal to the desired state. Suppose, 
for example, we wish to make an even cat after detecting 
to = 2, but three photons actually arrive at the detector. 
If the detector registers only two photons, then we have 
produced an odd cat, but we falsely believe we have an 
even cat. Other researchers have analyzed photon de- 
tector inefficiency in photon subtraction for the case of 
subtracting only one photon or using photon detectors 
that cannot discriminate 1, 2, ... photons [H, [H, l60l |. 
Here we will consider true photon counting. Let the de- 
tector's efficiency be represented by rj. The probability 
to register to photons (meaning that the detector reports 
seeing to photons, when more than m may actually arrive 
at the detector) is given by Eq. (fl"4|) . The state produced 
when the detector registers to photons is now given by 
the density operator 



' v ' n=m 

(56) 

The fidelity of p v (m) depends on the probability dis- 
tribution given by P(n). In fact given that to pho- 
tons are registered, the fidelity will be greatest when 
P(m + 1) -c P(to). When this condition holds, we can be 
confident that the detector's registering to photons is not 
caused by the arrival of more than to photons and the 
detector's inefficiency. This condition can be achieved by 
using less squeezing and a beam splitter with larger T. 
Of course, this will result in a decrease in the probability 
to register to photons. By choosing the best combina- 
tion of A and T we can obtain a fidelity that is equal 
to that achievable with perfect detectors. Unfortunately, 
with the optimal combination of A and T, P(to) is very 
small. Let us continue our examination of the a = 2 
case, but now the total efficiency of the photon coun- 
ters and other optical elements in the photon counting 
mode is 77 = 0.9. We would like to understand the trade- 
offs between fidelity and probability for various photon 
detection goals. In figures [T2l HSl and [14] we plot the 
fidelity, P^(to), and the beam splitter transitivity T as 
functions of A for m = 2, 4, and 6. In each case T has 
been optimized to give the highest fidelity. In each case, 
the maximum fidelity is achieved in the a region of low 
P,,(to). The maximum fidelity in the presence of loss is 
equal to the maximum fidelity achievable when perfect 
photon counters are used. 

We now expand our analysis to include two additional 
error sources: impurity in the squeezed state and "dark 
counts" in the photon counter. Ideally, the squeezed 
state would be a pure state of a mode matched to 
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FIG. 13: Probability to detect m = 2 (solid curve), m — 4 
(small dashes), and m = 6 (long dashes) as a function of A 
when detectors with efficiency r/ = 0.9 are used. For each com- 
bination of A and m, the beam splitter transitivity has been 
adjusted to give the maximum fidelity, as shown in Fig. 1121 




FIG. 14: Beam splitter transitivity T that gives the largest 
maximum fidelity in Fig. [12] as a function of A. The solid 
curve shows the m — 2 case, the small dashed curve shows 
the m = 4 case, and the large dashed curve shows the m — 6 
case. 



the local oscillator. However achieving this goal is ex- 
tremely difficult, especially for pulsed squeezing experi- 
ments [6l|, [62], |63j, [6J] . Instead, the squeezed state in the 
mode matched to the local oscillator is mixed, and other 
unwanted photons are traveling down the same beam 
path in other modes (for example, the wave packet car- 
rying the extra photons has a different temporal profile) . 
These unwanted photons may be registered by the pho- 
ton counter, falsely signaling the creation of a cat state. 
For these reasons we expect that these unwanted pho- 
tons will behave much like "dark counts." The homo- 
dyne detection system will be immune to contamination 
by photons in non-matched modes, because homodyne 
detection is sensitive only to the mode exactly matched 
to the local oscillator. 

We use the model in Fig. [15] to incorporate squeezed 
state impurity, photon counter inefficiency, and "dark 
counts." The squeezing S±(r) creates a pure squeezed 
state in mode 1. Mode 1 then passes through a beam 
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FIG. 15: Model of photon subtraction scheme including im- 
purity in the initial squeezed state, inefficiency of the photon 
counter, and "dark counts." The beam splitters with trans- 
missivities v and rj model loss from the squeezed state and 
inefficiency in the photon counter. Modes 2 and 4 are lost 
to the environment, mode 5 contains "dark counts", and the 
output of mode 1 contains pdijn), which should be similar to 
a cat state. 



splitter with transmissivity v, and 1 — v of the light is 
lost to the environment. In any real squeezing experi- 
ment, the ideal squeezing and decoherence takes place 
simultaneously inside the down-conversion crystal. This 
model of ideal squeezing followed by a beam splitter can 
create any state whose x— and p— quadratures are Gaus- 
sian and centered at the origin, provided that the vari- 
ance of one of the quadratures is equal to or less than the 
vacuum's variance. The observed quadratures' variances 
are related to the squeezing r and v by 



l-v 



2 ' 



and 



-v e 2r 

h v . 

2 2 



(57) 



(58) 



We use the conventions that the vacuum's variance is 1/2 
and h = 1. 

The impure squeezed state in mode 1 then encounters 
the beam splitter with transmissivity T, which diverts 
light into mode 3. To model the photon counter, mode 3 
passes through the beam splitter, whose transmissivity rj 
is equal to the real photon counter's efficiency. Mode 4 
is lost to the environment. In addition to mode 3, mode 
5 also enters the photon counter. The probability that 
mode 5 contains q photons is given by the Poisson distri- 
bution with mean d, given in Eq. (|15p . The actual "dark 
count" probability distribution may be different from the 
Poisson distribution, especially when many of the "dark 
counts" are actually caused by unwanted photons gener- 
ated in the squeezing. 



15 



We calculate exactly the full state (that is, the x- 
quadrature wave functions) of the five modes symboli- 
cally with the Mathematica computer algebra software. 
The output of mode 1, conditioned on m photons be- 
ing registered by the photon counter, is the mixed state 
Pd(m). We calculate the fidelity between pd(m) for typi- 
cal values of v, rj, and d. We would like to optimize the 
fidelity by adjusting r, T, and m. It is somewhat unre- 
alistic to treat r, is, and d as independently adjustable 
parameters, because they are all governed by dynamics 
occurring inside the down-conversion crystal. However, 
this model should still be useful for making some general 
recommendations and predictions for experiments. 

We consider two representative imagined experiments. 
For the first, whose goal is to make an odd cat with 
a = 1.25 by measuring m = 1, we have v = 0.85, 
7] = 0.8, d = 0.0005. The maximum achievable fidelity 
with these parameters is F = 0.727 using T = 0.922 
and r = —0.514. After the pure squeezed state passes 
through the beam splitter this corresponds to squeezed 
quadrature variances of 4.2 dB above the vacuum and - 
3.5 dB below. The detector will register one photon (or 
one "dark count") with probability 0.0143. 

For the second experiment, the goal is to make an even 
cat with a — 2 by measuring m = 2. For this more 
ambitious effort, v = 0.95, i] = 0.9, and d = 0.0002. 
The best fidelity is now F = 0.737 using T = 0.982 and 
r = —0.722. This corresponds to the observable squeezed 
quadrature variances of 6.11 dB above the vacuum and - 
5.62 dB below. The detector will register two photons (or 
one photon and one "dark count" or two "dark counts" ) 
with probability 0.00021. 

In the presence of significant levels of "dark counts," 
the strategy of reducing \r\ and increasing T to compen- 
sate for detector inefficiency is no longer effective. There 
is now a maximum m such that conditioning cat produc- 
tion on larger numbers of photons creates a lower fidelity 
cat state. High purity of the initial squeezed state and 
matching of the squeezed state's spatial and temporal 
mode shape to the local oscillator is critical for making 
high fidelity cat states. Any photons produced in the 
squeezing process that occupy modes other than that 
matched to the local oscillator will pollute the photon 
counter, degrading the cat fidelity. It may be necessary 
to engineer the pump mode shape and the squeezing crys- 
tal to reduce the squeezing of other modes using methods 
similar to those described in [65l 1661 167| • 

In recent years some experiments have used photon 
subtraction to make small cat states, which are some- 
times called "Schrodinger kittens" . All of these exper- 
iments used standard APDs to subtract the photons. 
These detectors are unable to discriminate 1,2,3,... pho- 
tons, but the squeezing level is low, so the probability of 
more than one photon being present is very small. 

The first of these experiments was performed by 
Wenger, Tualle-Brouri, and Grangier [lj|. They used 
squeezed light of r = —0.430, from which they extracted 
a single photon using a beam splitter with T = 0.885 



to create an odd cat state. They used a laser pulsed 
at 790 kHz and made 30,000 cat states in a three hour 
experiment. This corresponds to a Pd(l) ~ 3.5 x 10~ 6 
or making ~ 2.77 cats/second. They used tomography 
to reconstruct the Wigner function of this state, which 
clearly shows non-Gaussian characteristics. However, the 
reconstructed Wigner function did not contain any neg- 
ative values (a feature common to all cat states). This is 
primarily because their experiment suffered from a high 
"dark count" rate, low photon detector efficiency, and 
impurity in the squeezed state. A second pulsed experi- 
ment was performed by Orjoumtsev, Tualle-Brouri, Lau- 
rat, and Grangier [HI ]. This improved on the previous 
experiment by having more careful filtering of the mode 
measured by the APD and higher purity in the initial 
squeezed state. They produced cat states with a = 0.89, 
0.84, and 0.78 having fidelities of 0.70, 0.64, and 0.63. 

A continuous wave experiment was performed by 
Neergaard- Nielsen and co-authors [l||. Because light 
flows through their system continuously, but the photon 
counter gives discrete clicks, they recorded 1000 homo- 
dyne measurements lasting 2 (is before and after each 
APD click. The series of homodyne measurements was 
then turned into a single quadrature value by weighting 
each data point according to a temporal mode function 
optimized to give the most negative Wigner function of 
the reconstructed state. In this way they were able to 
effectively extract a single pulsed mode from the con- 
tinuous wave. They produced cat states with a = 1.05 
and 1.3 with fidelities of 0.53. A similar continuous wave 
photon subtraction experiment was performed by Wakui, 
Takahashi, Furusawa, and Sasaki [l7|]. They showed 
clear evidence of nonclassical states with negative Wigner 
functions, but they did not report cat state fidelities. 

Progress in making higher fidelity and larger amplitude 
cats by photon subtraction can be achieved by carefully 
tuning the beam splitter transmissivity and the level of 
squeezing. Until now, all experiments have subtracted 
only one photon from the squeezed state. Using pho- 
ton counters that can discriminate 0, 1, 2, ... photons 
would allow one to reject multi-photon events that de- 
grade the fidelity of the cats produced in previous ex- 
periments. Higher fidelities and amplitudes can also be 
achieved by subtracting more than one photon. Sub- 
tracting two or more photons allows one to use a less 
squeezed state, which can be produced with higher pu- 
rity. However, the multimode nature of pulsed squeezed 
states may have a greater deleterious effect when more 
photons are subtracted. 



VI. MAKING KITTENS 

A small odd cat state is well approximated by a 
squeezed single photon [ljj . When the squeezing op- 
erator S(r) = e^ a _a+ ) is applied to a single photon, 
the resulting state can be expanded in terms of photon 
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FIG. 16: Fidelity between an odd cat state of amplitude a 
and a squeezed single photon as a function of a. 




FIG. 17: Squeezing r required to maximize the fidelity be- 
tween a squeezed single photon and an odd cat of amplitude 



number states as 



n=0 



(tanhr) J V(2?i + 1)! 
(coshr)^2"n! 



\2n+l), (59) 



where r is the squeezing parameter. The fidelity of this 
state with an odd cat state is given by 



F(r,a) = |<*_(«)|5(r)|l)| 2 

2a 2 exp[ev 2 (tanh r — 1) 



(coshr) 3 (l — exp[— 2a 2 ]) 



(60) 



The maximized fidelity as a function of a is shown in 
Fig. I16[ll] . and the r that provides the maximum fidelity 
is shown in Fig.[T7|68|. The odd cat is well approximated 
by a squeezed single photon for small amplitudes a and 
small amounts of squeezing r. A single photon squeezed 
by r — 0.31 would have a fidelity of 0.997 with an odd 
cat of amplitude a = 1. 

There are numerous technologies for producing single 
photons. Some of the most promising include heralded 
photons from nondegenerate down-conversion, optically 
and electrically excited quantum dots, and atoms and 



FIG. 18: Fidelity of the odd cat state with a = 1/2 with 
squeezed single photon state made from a photon source 
whose inefficiency is p. 



ions in optical cavities. We refer the reader to 69] for pa- 
pers on single photon sources. Unfortunately, pure single 
photon states are not produced using current technology. 
The single photon is always in a mixture with the vacuum 

as 



p|0)(0| + (l- P )|l)(l| 



(61) 



where p is the inefficiency of the photon production. For 
most schemes, the probability to produce two or more 
photons is usually much smaller than p, so we ignore that 
possibility here. Therefore, the squeezed single photon 
state is also mixed with a squeezed vacuum 

pS(r)\0)(0\SHr) + (1 - p)S(r)\l) (l\&(r). (62) 

The squeezed vacuum is orthogonal to an odd cat state, 
so the fidelity of the state given by (|62|) and an odd cat 
state is (1 — p)F(r, a). We plot the fidelity of the state in 
Eq. (|62|) in Fig. [18] as a function of the photon production 
inefficiency p for a = 1/2. 

The success of this scheme to make kittens requires 
improvements in single photon production technologies. 
It also requires a practical method to coordinate the sin- 
gle photon production with squeezing. One could at- 
tempt such an experiment by continuously pumping the 
down-conversion crystal, but homodyne measurements 
are recorded only when a single photon is (or is supposed 
to be) present. 



VII. GROWING KITTENS INTO CATS 

Lund, Jeong, Ralph, and Kim described a method for 
producing larger amplitude cat states given a source of 
kittens [llj . This process uses only linear optical devices 
and photon detectors, and it is resilient to detector inef- 
ficiency. The scheme is depicted in Fig. [19J 

We begin with two small amplitude cat states \^^(a))i 
and \^ v ((3))2 in modes 1 and 2. The kittens meet in the 
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FIG. 19: Schematic of the nondeterministic amplification pro- 
cess to grow cat state amplitude. Two small cat states at 
modes 1 and 2 and a coherent state in mode 3 are manipu- 
lated with linear optics. A larger amplitude cat is produced 
in the output of mode 2 when photons are detected in both 
modes 1 and 3. 



first beam splitter, whose transmissivity is set to 

.,2 



T 



a 



(3 2 



(63) 



One output of that beam splitter is sent to a second beam 
splitter, which has transmissivity 1/2. Mode 3 contains 
the coherent state 



|7>3 = 



2af3 



^a 2 + P 2 



(64) 



The two beam splitters transform the two small cats and 
the coherent state as 



B 13 (l) B 12 (T)|^(a)) 1 |^(/3)) 2 | 7 ) 3 cx 
+e^|V2 7 > 1 



a 2 -[3 2 



+e^|0)i 



a 2 -0 2 



(65) 



where A = y/ a 2 + (3 2 . When modes 1 and 3 exit from the 
second beam splitter, they are sent to photon detectors. 
It is clear from Eq. (|65|) that the resulting state for mode 
2 when both detectors register photons is a large cat oc 
(I -A) 2 +e^f + ^\A) 2 ), whose amplitude A = ^a 2 + (3 2 
is larger than both a and (3, and the relative phase is the 
sum of the relative phases of the input small cats. 

The success probability (that both detectors detect at 
least one photon) for a single attempt of the process 
above is 
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FIG. 20: Success probabilities for a single attempt at the 
process to make a cat from two kittens, depicted in Fig. [19] 
We show probabilities for the input fields of two identical 
small odd cat states (solid curve), two identical small even 
cats (dashed curve), and even and odd small cats (dotted 
curve) . 



which we plot for different input combinations in 
Fig. |2"0"[ 1 1] ■ The probabilities depend on the type of small 
cats (odd or even) used. The probability P v ^(a,a) for 
two identical odd cats is always larger than 0.214, regard- 
less of the value of a. 

One could make the small cats necessary to produce a 
larger cat by this scheme, using any of the methods in 
this paper. 

The photodetector's inefficiency will not affect the fi- 
delity of the produced state. It reduces only the success 
probability. However, the detectors' dark counts may 
cause errors. Events in which the detector registers more 
photons than actually arrive at the detector can cause us 
to accept a state that we should discard. If the real num- 
ber of photons coming to either of the detectors is zero, 
we should discard the resulting state, but dark counts can 
make the detector register a number of photons different 
from zero and the resulting state would be accepted, de- 
creasing the fidelity of the produced state. On the other 
hand, the inefficiency of the detector decreases the prob- 
ability of generating the correct state. Here we provide 
new analysis of this scheme incorporating detector inef- 
ficiency and dark counts. 

Although inefficiency alone would not affect the fi- 
delity, to correctly model the effect of dark counts we 
must also include the photon detectors' inefficiency. We 
imagine a beam splitter in front of each detector with 
transmissivity rj. We suppose that l\ and I3 photons pass 
through the beam splitters before the detectors in modes 
1 and 3. This gives the state 



P v (h,h) 



ni—h 713-/3 

XT) h (l-ri) ni - h x r] h (l-riy i:, - h 
x|V>(rci,n3))(V>(ni,7i3)|, (67) 



(66) where ^(ni, 713)) is the normalized state given by the 
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FIG. 21: Fidelity of a cat with a = 2 produced from kittens 
with ai = y/2 as a function of the average number of dark 
counts in the detectors. The detectors have efficiencies rj — 
0.88 on the solid line and r\ = 0.8 on the dashed line. 



FIG. 22: Fidelity of a generated cat with a — 2 produced 
from kittens with at = \p2 as a function of the detectors' 
efficiency r\. The detectors have mean dark counts d — 10 -4 
on the solid curve and d = 10 -2 on the dashed curve. 



projection of Eq. (|65p into (rti , 713 1 , and P(ni,n 3 ) is the 
probability of measuring n\ and n 3 photons in modes 
1 and 3, respectively. To construct the density matrix 
of the output mode after the dark counts we must add 
the Pr/ihjh) with their correct probability weights. This 
gives us 

^ mi m.3 

p d (mi,m 3 ) = -5-7 r V V P v {h,h) x 

p<i(mi - h)vd(m 3 - h)pri(h,h),{6&) 

where P^(/i,Z 3 ), pd(m,i - Zi), and P d (mi,m 3 ) are given 
by Eqs. (fT4|) , (fl~5|) . and (fl~7)) , respectively. j0d(mi,m 3 ) is 
the state produced when the detectors register mi and 
m 3 photons and dark counts. The scheme is believed to 
have succeeded when both detectors register one or more 
photons, so the full state produced by this scheme is 

Pacccpt — v~>°o n / \ ' V""J 

L TOl =iL m ,=iPd(mi.n>8) 

The fidelity is F = (\E , ±(a)|/3acce P t|^'±(a))- Let us con- 
sider the generation of a larger cat by this amplification 
scheme, given that we have two small odd cats as input. 
If we use a photon counter whose efficiency is rj = 0.8 and 
an average number of dark counts per detection event of 
d — 4x 10~ 4 , we can generate an a — 2 cat with fidelity 
F = 0.9994. If we use a TESPC, described in section [fl 
whose efficiency is r\ = 0.88 and the average number of 
dark counts per detection event is d = 10~ 8 , we can ob- 
tain an a = 2 cat with a fidelity of F = 0.999999986. To 
obtain such an average number of dark counts we need to 
carefully filter out blackbody radiation, shield the exper- 
iment from ambient light, and restrict detection events 
to short time windows. 

Figs. and fJU show the fidelity as a function of the 
average number of dark counts and the efficiency, for gen- 
erating an a — 2 cat state. 



This scheme has been analyzed in more detail by 
Jeong, Lund, and Ralph in [70]. An important property 
they discussed is that, when using impure input kittens, 
this process can result in a purer output cat state. That 
is, the output state can have a higher fidelity with a true 
cat than the fidelity of the input states. Of course this is 
not generally true for any input state, but it is very effec- 
tive when the input states are squeezed single photon kit- 
tens produced from inefficient single photon sources, as 
we discussed in the previous section. This is because the 
cases in which a squeezed vacuum rather than a squeezed 
photon is input are more likely to be rejected because one 
of the photon detectors does not click. 

According to the calculations in [70] , if the probability 
that a single photon is not present is p — 0.4 and the in- 
put squeezed single photons approximate cat states with 
amplitude a = 1/2, the fidelity of the initial small cat is 
F = 0.60, but the fidelity of the output state after one 
iteration is F = 0.89. If p = 0.25 {p = 0.05), the fidelity 
of the initial cat state is F = 0.750 (F = 0.950) but be- 
comes F = 0.941 (F — 0.990) after one iteration of the 
cat growth procedure. 

Jeong and co-authors also report that when one gen- 
erates a large cat out of perfect squeezed single photon 
states, the fidelity maximizes for a particular number of 
iterations. Starting with smaller ai states and perform- 
ing more iterations to achieve a = 2 would produce a 
lower fidelity. For example, four iterations starting from 
the initial amplitude a, = 1/2 are required to gain the 
maximum fidelity F — 0.995 for a = 2. High fidelity, 
F > 0.99, can be obtained up to a = 2.5. 

A slight adaptation of the scheme in Fig. [T9l was pro- 
posed by Rhode and Lund in [7l|. They suggested using 
one even cat and one odd cat as input states. The even 
and odd cats interfere at a beam splitter. One of the out- 
put ports of the beam splitter is measured by a photon 
detector. A larger amplitude odd cat state will appear 
in the other output port when the detector detects zero 
photons. Their paper also contains some analysis of the 
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FIG. 23: Schematic diagram of the process to generate a cat 
state using small Kerr nonlinearity, a beam splitter, and a 
homodyne detection of the i-quadrature. 



performance of this scheme in the presence of mode mis- 
match and loss. 

Given a source of kittens, this appears to be an ex- 
cellent scheme for producing larger amplitude cat states. 
It can be implemented with current technology. Its re- 
silience to photon detector inefficiency and dark counts 
and its ability to purify its input are particularly attrac- 
tive characteristics. 



VIII. SMALL KERR EFFECT 

In 2004, Jeong and co-authors [l|| proposed a proba- 
bilistic scheme to generate cat states using a small Kerr 
effect, a beam splitter, and a homodyne measurement. 
This scheme uses a self-phase modulating Kerr medium 
like that discussed in Section [D] but the Kerr Hamilto- 
nian is applied for a much shorter time, so photon loss 
should be a less significant factor. A diagram of the 
scheme appears in Fig. |2"3"1 

We begin the process by subjecting an initial coherent 
state in mode 1 to the Kerr Hamiltonian given by 
Eq. (|2"4")l . Under the influence of this Hamiltonian, an 
initial coherent state will evolve to the following state 
after time t: 



(70) 



where 4> n = %in 2 . The resulting state in Eq. ([70)) is an 
example of a generalized coherent state introduced by 
Titulaer and Glauber [lH and discussed by Bialynicka- 
Birula [73]. Generalized coherent states can always be 
represented as superpositions of coherent states, with the 



same modulus but different phases. In the specific case of 
generalized coherent states for which 0„ is periodic in n, 
they can be represented as a discrete superposition of N 
coherent states. For example, when the interaction time 
Xt is ir/N with a positive integer N, the initial coherent 
state |a)i evolves to [74[ 
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where we neglected the normalization. C n ^ obeys 



JV 



2 i n it /N\h _ e -iizk 2 /N 



(72) 



n=l 



Using the method of Gantsog and Tanas [75j , this can be 
solved to determine that 



JV-l 
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(73) 



k=0 



The state in (I71[) is now combined with the vacuum 
state |0)2 in mode 2 on a 50-50 beam splitter. After 
passing through the beam splitter the system is in the 
state 
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We next perform a homodyne measurement of the x- 
quadrature of mode 2. The measurement reduces the 
state of Eq. (J74J) to 
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C n ,N{x) = N x C n<N 2 (x 
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N x is a normalization factor, and |a;) is the observed x- 
quadrature eigenstate. To obtain the desired cat state, 
the measurement result is selected in such a way that 
the coefficients \Cn/2,n{x)\ and \Cn.n(x)\ have the same 
nonzero value and all the other coefficients are zero. For 
example, if x = is measured and N = Ah, where k is an 
integer, the coefficients will be the largest for n = N/A 
and n = 3N/A. The coefficients will become smaller as n 
gets far from those two points. 

The fidelity between the state \4>n(x)} and a perfect 
cat with appropriate amplitude is given by 



F(ai, N, x) — max 



;* (a I /V2)|V'jv(x)) 



(77) 



20 




10 




FIG. 24: Fidelity F of the cat state produced by a weak Kerr 
effect when an x = measurement is obtained versus the 
generated cat amplitude ia. The solid curve shows the case 
when N — 20 (one tenth of the interaction required by using 
the Kerr effect directly), the small dashes show N = 12, and 
the long dashes show N = 8. 



Since the scalar product of two coherent states is 



e -|/3-a| 2 /2 e (a*/3- a /3*)/2 



(78) 



we can write the fidelity as [lj 



F = max 



exp 



N 



n=l 



e" v exp 



4( 



1 - e 



(79) 



For an initial coherent state a, and an interaction time 
Xt = 7r/20 (N — 10, one tenth of the interaction required 
to make a cat state using only the Kerr medium), the 
fidelity between the state that leaves the Kerr medium 
and an ideal cat is F s» 0.1 (for ati — 20). After the 
homodyne measurement, the reduced state is similar to 
a cat state with amplitude a = ion /y/2. The maximum 
fidelity is achieved when x — is observed. The output 
state is most similar to an even cat {<j> = 0) state. Figure 
[Ml shows the fidelity of a generated cat when x = is 
measured as a function of the final amplitude a. A good 
fidelity is obtained only for a final amplitude a larger 
than « 5 when N = 20. 

The highest fidelity is achieved for the measurement 
result x — 0. We show a plot of the fidelity as a function 
of x in Fig. [55] The fidelity shown in the plot is the 
maximum fidelity for all possible values of <f>. The phase 
of the cat state produced oscillates as a function of x , as 
shown in Fig. [^Hl .76]. Each cat's phase is random but 
can be inferred from the measurement of x. Although 
the random phase is inconvenient, because it is known, 
one may correct or account for it. 

To ensure high fidelity cats, one would choose some 
5 such that only measurement cases with |a;| < S are 



FIG. 25: Fidelity of the generated cat state against the mea- 
surement outcome x for N = 20. The solid curve shows gen- 
erated cat amplitude a = 20i, and the dashed curve shows 
a = lOi. The fidelity has been maximized over the phase (j> 
of the generated cat state. 
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FIG. 26: Cat state phase <f> that maximizes the fidelity shown 
in Fig. [25] The solid curve shows generated cat amplitude 
a — 20 i, and the dashed curve shows a — Wi. 



accepted. The success probability to get an acceptably 
small x is fijf 



,N,S) 



dxTru 0^)1212(^1 <%> \x)n(x\ 



J ° n,m=l 



x exp 



a i /-1 _ 2i(m-n)n/N\ 

' 2 1 1 



, (80) 



where \S\is the range of accepted ^-quadrature measure- 
ments. This "success probability" is the probability to 
obtain a cat of any phase </>. The probability to make a 
cat with a particular phase is much smaller. 

One can make cat states with a = 20i using N = 20 
and accepting all measurement results with |x| < 3.75. 
All of these states would have fidelities greater than 
0.99997 with perfect cats with correctly chosen phases. 
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FIG. 27: Diagram of a scheme to produce cat states using a 
cross-phase modulating Kerr medium. The dashed lines are 
beam splitters with transmissivity equal to 1/2, and the solid 
bold lines are mirrors. The triangles are photon detectors. 

The probability to obtain \x\ < 3.75 is 0.052. To make 
cats with a = 10% and that have fidelities greater than 
0.9991, one should accept |x| < 1.06, which will happen 
with probability 0.045. When making cat states with an 
even shorter interaction time (larger N) , the fidelities are 
smaller and the success probability is lower. 

A related scheme for using the Kerr effect was proposed 
by Gerry This scheme uses a cross-phase modulat- 
ing type of Kerr medium, a single photon source, and 
photon detectors. We picture the scheme in Fig. [27] 

A single photon is inserted into a Mach-Zehnder inter- 
ferometer, occupying modes 2 and 3. The coherent state 
| a) is in mode 1. After the photon enters the interferom- 
eter, the three modes are in the state 

^(M,0)i 23 + k0,l) 123 ). (81) 

One arm of the interferometer (mode 2) and the coherent 
state in mode 1 pass through a Kerr medium, governed 
by the Hamiltonian xfi 1 n 2 . After time t, this interaction 
evolves the system into the state 

-^=(|ae-^,l,0) 123 + |a,0,l) 123 ), (82) 

where cf> — \t- The phase of the coherent state is 
shifted when the photon accompanies it through the Kerr 
medium. After modes 2 and 3 pass through the second 
beam splitter and exit the interferometer, the system is 
in the state 

i(|ae-^,l,0) 123 + |ae-^,0,l) 123 

-|a,l,0)i 23 + |a,0,l)i 23 ). (83) 

When detector A sees (1) photons and detector B sees 
1 (0) photon, the resulting state of mode 1 is 

4=(|ae- # )i±Hi). (84) 



When 4> — Xt = n this produces an even or odd cat 
with equal probability, depending on which of the two 
detectors registers the photon. 

Like the method discussed in Section HH this requires 
a Kerr medium with a large ratio of nonlinear strength 
to loss (x/l) to make high fidelity cats. However, in [77| 
Jeong proposed an adaptation of this scheme designed 
to overcome the problem of loss in the Kerr medium. 
Suppose one desires to make the cat state of amplitude 
a. Rather than inserting the coherent state \a) into the 
network in Fig. [27l one should insert a coherent state 
with a much larger amplitude (3 » a. One then 
applies the Kerr effect for a much shorter time, so that 
\f3e %xt — j3\ = 2a. The displacement operator can shift 
the state \(3e~ txt ) ± \/3), so that it is symmetric about the 
origin of phase space. This results in the desired state 
\a)±\-a). 

Jeong also analyzed the behavior when the Kerr 
medium absorbs photons. If we increase /3, the coherence 
of the desired output state, of fixed amplitude a, also in- 
creases. This is possible because one can subject the light 
to the lossy Kerr medium for a shorter time. According 
to Jeong's estimates one could make a high fidelity cat 
with a = 3 using /3 ~ 30, 000 through an optical fiber of 
190 m. While Jeong's calculations show that this scheme 
is robust against loss in the Kerr medium, it requires dis- 
placing a superposition of two coherent states of many 
photons to the origin of phase space. Small displace- 
ments can be achieved with high fidelity by passing the 
state through a beam splitter with transmissivity T — ► 1, 
while a high amplitude coherent state passes through the 
other input of the beam splitter. It is not clear how well 
this displacement scheme will work for states contain- 
ing ~ 30,000 2 photons. Any phase instability or mode 
mismatching in the displacement operation may lead to 
significant errors. 



IX. CONCLUSIONS 

We have reviewed and analyzed several methods to 
make cat states. Aside from making the cats directly 
using a strong, low loss Kerr effect, all of these methods 
involve preparation of the light using a lesser nonlinearity 
and postselecting on some measurement. We attempt to 
compare some of the characteristics of the schemes in 
Table El 

We cannot make a definitive recommendation for the 
"best" cat making scheme for all circumstances. This is 
because the schemes were formulated with different ob- 
jectives, such as desired cat state amplitude, fidelity, and 
success probability; also the schemes each require signifi- 
cantly different resources and technologies to implement. 
Nevertheless we can make a few comments on the relative 
utility of the schemes. 

The photon subtraction scheme seems to be the sim- 
plest to implement for a laboratory demonstration of cre- 
ating small amplitude (a < 2) cats. In fact a few exper- 
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TABLE IV: Comparison of the Features of Cat Making Schemes 





Detectors 


Nonlinearity 


Cat Produced 


Direct Kerr effect 


None 


Strong, low loss Kerr effect 


F — 1 for zero loss 


Back-action evasion 


Efficient photon counting 


Squeezing twice 


High fidelity, low probabil- 
ity, low amplitude 


Photon subtraction 


Efficient photon counting 


Degenerate squeezing 


High fidelity, low probabil- 
ity, low amplitude 


Squeezed single photon None 


Single photon input, squeez- 
ing 


Very low amplitude 


Kittens into cats 


Photon detection low effi- 
ciency is OK 


None 


F — 1 for perfect input kit- 
tens 


Small Kerr effect 


Homodyne 


Brief or weak Kerr effect 


High amplitude, fidelity con- 
ditional on homodyne mea- 
surement 



iments have already succeeded in this task [lfj, llj} • Im- 
proving on these experiments is likely to require making 
squeezed states with greater purity (especially in pulsed 
experiments) and subtraction of more than one photon. 
Given the simplicity and fidelities achievable with the 
photon subtraction scheme, there seems to be little use 
in implementing any of the back-action evasion schemes. 
The back action evasion schemes provide only slightly 
higher fidelities but require two stages of squeezing and 
are significantly more complex to implement. 

If a large number of high fidelity small amplitude cats 
is required for use in a quantum computer, it may be 
advantageous to make the cats by squeezing single pho- 
tons. In this application it will be helpful to have the 
cats available on demand, which will require a high fi- 
delity triggered single photon source, such as those based 
on quantum dots. Quantum computation with cat states 
is likely to require fidelities significantly larger than 0.99, 
perhaps even 0.999 - 0.9999. The iterative scheme for 
growing kittens to cats may serve to both increase cat 
amplitude and fidelity. The combination of making kit- 
tens by squeezing single photons and then growing the 
kittens has the advantage that the resulting cat fidelity 
does not depend on the efficiency of photon detectors. 

High amplitude cats may be most useful for interferom- 
etry applications. In the future, making high amplitude 
cats might best be accomplished using the Kerr effect, 
either directly or after postselecting after homodyne de- 
tection as described in the "Small Kerr Effect" section. 
Making high fidelity cats directly through the Kerr ef- 
fect requires a loss to nonlinear strength ratio j/x ^5 !• 
Single mode fused silica fibers have j/x ~ 260, and we 
have been unable to find any material with a smaller 
ratio. Producing cats in this manner will likely require 
significant advances in technologies, such as EIT or pho- 
tonic crystals. Proposals for using a weaker Kerr effect 



followed by postselection may have less strict require- 
ments for 7/x, but they are effective only when used 
with large numbers of photons, so large that the cats will 
be extremely susceptible to decoherence due to photon 
absorption and other errors. 

Much previous research in Kerr materials has been de- 
voted to increasing the interaction strength while suffer- 
ing from much larger photon loss. We strongly encour- 
age research devoted to making low loss Kerr materials 
or interactions. Low loss Kerr media will have many uses 
in addition to making cat states such as non-demolition 
photon counting and a controlled-phase quantum logic 
gate for single photon qubits. 

Mode mismatch is an error source that can affect any 
of the methods to make cat states. While mode mismatch 
sometimes behaves similar to photon loss, its effects can 
be more subtle. This is an issue that deserves more anal- 
ysis than we have provided here. 
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